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Abstract 

The consistency of the constraint with the evolution equations for spatially inhomogeneous and 
irrotational silent (SIIS) models of Petrov type I, demands that the former are preserved along 
the timelike congruence represented by the velocity of the dust fluid, leading to new non-trivial 
constraints. This fact has been used to conjecture that the resulting models correspond to the 
spatially homogeneous (SH) models of Bianchi type I, at least for the case where the cosmological 
constant vanish. By exploiting the full set of the constraint equations as expressed in the 1+3 
covariant formalism and using elements from the theory of the spacelike congruences, we provide a 
direct and simple proof of this conjecture for vacuum and dust fluid models, which shows that the 
Szekeres family of solutions represents the most general class of SIIS models. The suggested procedure 
also shows that, the uniqueness of the SIIS of the Petrov type D is not, in general, affected by the 
presence of a non-zero pressure fluid. Therefore, in order to allow a broader class of Petrov type 
I solutions apart from the SH models of Bianchi type I, one should consider more general "silent" 
configurations by relaxing the vanishing of the vorticity and the magnetic part of the Weyl tensor 
but maintaining their "silence" properties i.e. the vanishing of the curls of E a b,H a b and the pressure 
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1 Introduction 



The discovery of the Szekeres family of solutions of the Einstein's Field Equations (EFE), initiated a 
continuous research effort to understand the richness of the structure of this interesting class of inhomo- 
geneous models [l]-[6]. The main reason for this interest was that, due to the inhomogeneity condition, 
these solutions can in principle describe cither expanding or collapsing matter distributions leading for 
example to a deeper insight into the formation mechanism of the local structure inhomogeneities due to 
the amplification of small density fluctuations in the early universe [7]. 

Nevertheless, it emerged that Szekeres models were "hiding" a sufficient number of crucial, in some 
sense, properties in their intrinsic structure. At a geometrical level, it was proved immediately after 
their discovery that they admit no isometries i.e. there are no non-trivial Killing Vector Fields (KVFs) 
admitted by the spacetime manifold. Therefore one could reasonably expect that the absence of KVFs 
should make Szekeres models quite general within the class of inhomogeneous models. However, the 
well known kinematical properties of the vanishing of the acceleration and the vorticity of the fluid 
velocity (or, in the case of the vacuum subfamily of models, the timelikc vector field normal to the spatial 
hypcrsurfaccs t =const.), indicate that these models are more special than expected. In addition, Barnes 
and Rowlingson have shown that, apart from they being algebraically special (namely of Petrov type D), 
the magnetic part H ab of the Weyl tensor vanishes [8]. Consequently Szekeres models belong to a wider 
family of models satisfying H ab = which can then be of the Petrov type I, D or O. This family contains, 
as special cases, many well-known and important solutions like the standard Friedmann-Lemaitre (FL) 
cosmological model, the Spatially Homogeneous (SH) Bianchi type I models, as well as other cosmological 
or astrophysical solutions of the EFE with sound physical interest [9] . 

Furthermore this result opened the possibility that these models may represent the exact solution 
analogues of the inhomogeneous first order scalar perturbations of the FL model, since it has been shown 
that the magnetic part of the Weyl tensor does not contain scalar contributions at the linear regime 
[10]. Although it has been proved that this is not true at second (or higher) order perturbations, the 
undertaken analysis for models satisfying the conditions: 

p = = q a , 7r ab = 0, u a = 0, uj ab = 0, H ab = ^e ac e fC efbd u c u d = 0. (1.1) 

has uncovered a striking (and, from a mathematical point of view, very often desirable) feature of the 
complete set of the EFE: the decoupling of the spatial divergence and curl equations (or equivalently 
the constraints i.e. the fully projected derivatives normal to the timelike vector field u a ) from the 
evolution equations of the kinematical and dynamical variables [11]-[14]. As a result, there is no physical 
mechanism for gravity wave propagation {H ab = 0) or the presence of sound waves (j> = 0) between any 
nearby timelikc curves, thus leading to the name silent models. Furthermore, due to the above mentioned 
decoupling, the associated equations for the physical variables are reduced to a set of first-order ordinary 
differential equations (odes) which has been studied using methods from the theory of dynamical systems 
with a view to analyzing the asymptotic behaviour of cither expanding or collapsing collisionlcss matter 
configurations [15], [16]. 

However, one must complement the above discussion by examining whether a suitable set of initial 
data evolves consistently, which will ensure that Spatially Inhomogeneous Irrotational Silent (SIIS) models 
retain their generality within the class of inhomogeneous solutions of the EFEs. This is equivalent to 
demanding that the spatial divergence and curl equations (encoded in the set of the initial data), are 
consistent with the evolution equations hence, they are preserved identically along the timelike congruence 
u a without imposing new geometrical, kinematical or dynamical restrictions. Although it has been stated 
that the constraints are generally consistent [17], [18] however, it has been demonstrated that this is not 
the case [19]. In particular, the consistency analysis for the SIIS of the Petrov type I has shown that the 
constraints, coming from the silent conditions (1.1), are not trivial and after repeated propagation along 
u a , we expect to obtain an infinite chain of non-trivial algebraic relations which allows one to conjecture 
that there are no SIIS models of the Petrov type I [20, 21, 22, 23] 

Clearly, if we assume that the conjecture is true, the Szekeres family of solutions keeps a privileged 



2 



position within the SIIS models due to its geometrical and dynamical characteristics and can be exploited 
in order to clarify some open issues in spatially inhomogencous models. For example, the consistency 
of the Petrov type D models together with their non-symmetry property (i.e. the non-existence of 
isometries), indicate that in the case of the Petrov type I models the new constraints must be related 
with the existence of some type of symmetry. This expectation can be justified by the presence of the non- 
trivial algebraic (consistency) relations in the Petrov type I models which necessary lead to a reduction of 
the allowed models, similarly to the effect by the existence of a geometric symmetry. On the other hand, 
the qualitative study of the SIIS models shows that, if the Szekeres family of solutions is unique, we can 
restrict the dynamical systems analysis only to the Petrov type D invariant set [25]. Among others, this 
provides a way to reveal the geometric nature and to describe the dynamical role of the equilibrium points 
of the associated state space and examine to what extent the structure of the latter has similarities with 
the corresponding state space of the SH models or the inhomogeneous models with one spatial degree of 
freedom [25]- [39]. 

The above discussion initiated a study towards either proving the conjecture or, at least, presenting 
results that enforce its validity. For example, it has been shown the uniqueness of the vacuum Szekeres 
sub-family within the SIIS models and the non existence of Petrov type I dust models admitting a G3 
group of isometries, provided in both cases that the cosmological constant A vanishes [23, 24], [40]. On 
the other hand, the presence of a cosmological constant (i.e. the A— term generalization of the spatially 
inhomogeneous and irrotational models studied in [41]) cancels, in some way, the dynamical effects of 
the Petrov type I conditions, therefore under certain circumstances the conjecture would not hold and a 
broader family of models exists [42]. 

The purpose of this article is to revisit the consistency problem in vacuum and pressurclcss SIIS models 
and provide an alternative method to resolve it, using the time propagation of the constraints only up 
to third order and avoiding the appearance of complicated algebraic relations between the dynamical 
variables. The suggested approach is used to provide a conclusive answer to the consistency conjecture in 
the affirmative, in a transparent and fully covariant manner. Our method heavily relies on the use of the 
theory of spacclike congruences in conjunction with the 1+3 covariant formalism as applied in relativistic 
cosmology and astrophysics [43]- [47]. The paper is organized as follows: Section 2 is devoted to present 
the basic elements of the theory of timelike and spacelike congruences and subsequently the 1+3 and 
1+1+2 covariant description of the SIIS models, in order to express the resultant constraints in terms of 
the irreducible kincmatical quantities of the spacelike congruences. The main result of the paper is given 
in Section 3 where we prove the uniqueness of the Petrov type D family within the class of SIIS models. 
This is achieved by showing that models of the Petrov type I always admit a three dimensional Abelian 
group of isometrics with spacelike orbits, and thus reduce to the subclass of SH models of the Bianchi 
type I [48]. We should note that the above uniqueness result is still valid when we allow the presence of 
a non-zero pressure fluid configuration or a non- vanishing cosmological constant, provided in both cases 
that the extrinsic curvature of the spacelike hypersurfaces is not degenerated or, more generally, the shear 
eigenvalues arc not restricted to be proportional to the expansion of the timelikc congruence. Finally in 
Section 4 we draw our conclusions and discuss the ways of how the approach presented in the present 
article can be extended to check the broadness of more general silent models. 

Throughout this paper, the following conventions have been used: the pair (M., g) denotes the space- 
time manifold endowed with a Lorentzian metric of signature (—,+,+,+), spacetimc indices are denoted 
by lower case Latin letters a, b, ... = 0, 1, 2, 3, spacelike eigenvalue indices are denoted by lower case Greek 
letters a, (3, ... = 1, 2, 3 and we have used geometrised units such that 8irG/c 2 = 1 = c. 

2 Covariant description of SIIS models 

2.1 Elements from the theory of timelike congruences 

The 1+3 covariant formalism [47] starts by introducing a unit timelike vector field u a (u a u a = —1) which 
is identified with the average velocity of matter in the Universe, thus representing the congruence of 
worldlincs of the so-called fundamental (preferred) observers. The existence of the timclike vector field 
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u a generates, at any spacetime event, a unique splitting of all the geometrical, kincmatical and dynamical 
(tensorial) quantities into temporal and spatial parts using the projection tensor: 



Kb = g a b + u a u b , h a c h ch = h ab , h a b h b a = 3, h ab u b = 0. (2.1) 

Essentially, h ab represents the intrinsic metric of the 3-dimcnsional rest spaces of the observers u a , even 
if these 3-spaces do not constitute, in general, an integrablc submanifold S of M.. The external geometric 
structure of the submanifold S can be efficiently described in terms of the irreducible kincmatical parts, 
coming from the 1+3 decomposition of the first covariant derivatives of u a according to: 

g 

U a ;b = (?ab + ^h ab + UJ a b ~ U a U b (2.2) 

where 8 = u a ^ b h ab is the overall volume expansion (or contraction) rate, o~ ab = h a k h b l [u^-i) — is 
the shear tensor describing the rate of distortion of S in different directions (i.e. the change of its shape), 
^ab = h a h b ut^-n is the vorticity tensor and u a = u a[b u b is the four-acceleration vector field. Moreover, 
since u a corresponds to the average velocity of the matter fluid, the above kincmatical quantities have 
also a direct physical interpretation. Their dynamical behaviour can be studied using the Ricci identity: 

2«a:[6c] = Rdabc.U d (2.3) 

and taking into account the EFE. 

For the case under consideration the fluid pressure vanishes and the EFE are given by: 

Rab = ^ ( u aU b + hab) + ^9ab (2.4) 

where p is the energy density of the matter fluid and A is the cosmological constant. 

The complete 1+3 decomposition of (2.3) into temporal and spatial parts, leads to evolution and 
constraint equations of the kincmatical quantities. Because the timelike congruence u a is irrotational 
(u> a b = 0) and geodesic (u a = 0) and using the "silent" conditions (1.1), the Ricci identity reduces to two 
evolution equations for 9 and a ab and two constraint (divergence and curl) equations for the shear tensor 
[47]: 

h - la ' 2 <J ab a ab -\p + A (2.5) 
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0~ab = --Z®0- ab - a a C G cb - E a b + 77 (o dCT ) h ab (2.6) 



hacO kc , k = \h a H, k (2.7) 

^{a^mns^ ^ = (2 g) 

where 

E ab = C acbd u c u d (2.9) 
is the electric part of the Weyl conformal tensor and we have used the notation: 

K a ... = K a ... k u k (2.10) 

for the time derivative of any scalar or tensorial quantity. 

We should note that the shear-curl constraint (2.8) is the result of imposing the silent conditions (1.1). 
In order to obtain a closed set of equations one must also use the Bianchi identities 

Rab[cd-e] = 0. (2-11) 
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The associated 1+3 decomposition of (2.11) gives the energy and momentum conservation equations plus 
two evolution and two constraints equations for the electric and magnetic part of the Wcyl tensor which 
for the case of a SIIS model take the form: 



E ab = -^P°ab - 0E ab + 3E (a c a b)c - (E cd a cd ) h ab (2.12) 

p = -p9 (2.13) 

h ac E kc .^ k = l -h a k p, k (2.14) 

h K* e b)mns EnimUs = (2 15) 

e abcd E bk a c k u d = 0. (2.16) 



Similarly, the electric-curl constraint (2.15) follows from the silent conditions (1.1). The above closed sys- 
tem of equations completely describes the dynamics and can be used in order to determine the broadness 
of the family of SIIS models. This depends on the consistency of the constraints (2.8) and (2.15) with 
the set of evolution equations (2.5)-(2.6) and (2.12)-(2.13). Although equations (2.8) and (2.15) appear 
similar they have a completely different origin and one should expect that their consistency demands will 
be independent. Eventually this is indeed the case since it has been shown that the propagation of (2.8) 
does not generate further conditions and evolves consistently along the worldlines of the fundamental 
observers [21]. On the other hand propagation of (2.15) leads to a set of non-trivially satisfied equations, 
indicating that equation (2.15) represents a new constraint. It has been pointed out that the different 
character of (2.15) is a reminiscent of the implications of the fact that the vanishing of the evolution 
equation of a dynamical variable leads to a new integrability condition [22]. As we shall demonstrate 
in the next section, equation (2.15) is not, in general, consistent with the evolution equations (2.6) and 
(2.12)-(2.13) and represents the compact form of a set of "hidden" geometric restrictions. 

A direct consequence of equation (2.16) is that the shear a a b and the electric part E ab tensors commute 
which implies that they have a common eigenframe. As a result, if {x a , y a , z a } is a set of three mutually 
orthogonal and unit spacelike vector fields that constitute the spatial eigenframe of <j a b and E ab , we may 
write: 

E ab = Eix a x b + E 2 yaVb + E 3 z a z b (2-17) 
<r ab = aix a x b + a 2 y a yb + cr 3 z a z b (2-18) 
where E a ,a a are the associated eigenvalues satisfying the trace-free conditions: 

^£ a = ^V Q = 0. (2.19) 

a a 

Furthermore it has been shown that each of {x a , y a , z a } is hypersurface orthogonal [8] i.e. 

X[ a Xb;c] = V[aVb;c] = Z[ a Z b:c ] = 0. (2.20) 

which implies that x a — y a = z a = and local coordinates can be found such that the metric can be 
written: 

ds 2 = g ab dx a dx b = -dt 2 + e 2A dx 2 + e 2B dy 2 + e 2T dz 2 (2.21) 
where A(t, x, y, z), B(t, x, y, z), T(t, x, y, z) are smooth functions of their arguments. 
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2.2 Elements from the theory of spacelike congruences 

In many situations with clear geometrical or physical importance, apart from the existence of a preferred 
timclikc congruence, there also may exist a preferred spacclike direction x a (x a x a = 1) representing an 
intrinsic geometrical or physical feature of the corresponding model (for example, due to the existence 
of a spacclike KVF or the presence of a magnetic field). Consequently, a further 1+2 splitting of the 
3-dimensional space naturally arises leading to the concept of the 1+1+2 decomposition of the spacctime 
manifold. 

In complete analogy with the 1+3 decomposition, the starting point is to introduce the projection 
tensor: 

p a6 (x) = g ab + u a u b - x a x b = h a b - x a x b (2.22) 

P a c Pcb=Pab, Pa b p b a = 2, Pab u b = = Pab x b (2.23) 

which is identified with the associated metric of the 2-dimensional space X (the so-called screen space) 
normal to each vector field of the pair {u a ,x a } at any spacetime event. Then the geometric structure 
of X is studied by decomposing, into irreducible kincmatical parts, the first covariant derivatives of the 
spacelike vector field x a according to [43]- [46]: 



X a ;b = %,b (x) + yPafc + 7£ afc (x)+ X a X b - X a U b + p b C ± c U a + [2uj cb X C - iVfe(x)] U a 



where: 



Tab{*) 



PaPb 



X a; bP 



l 



:£xPkl 



%bP' 



a b 



K ab (x) = p a k p b l x [k .i] 



N a =p 



a (± k - U k ) 



(2.24) 

(2.25) 

(2.26) 

(2.27) 
(2.28) 



are the rate of the surface expansion, the rate of shear tensor, the rotation tensor and the Greenbcrg 
vector field of the spacclike congruence x a respectively and we have used the notation: 



K a .. 



■kX 



(2.29) 



for the directional derivative along the vector field x a of any scalar or tensorial quantity. 

The geometrical meaning of each of the defined kincmatical quantities of the spacelike congruence is 
similar to that of the corresponding quantities of u a , carrying information on the (overall or in different 
directions) distortion of X as measured by the fundamental observers u a . The new ingredient is the ap- 
pearance of the Greenberg vector iV Q which is of crucial importance in the theory of spacclike congruences 
and represents the "coupling" mechanism of the dynamical behaviour of the model in directions normal 
and parallel to the screen space X. For example, the equation N a = implies that the pair of vector 
fields {u a ,x a } generates a 2-dimcnsional intcgrablc submanifold of M. and the spacclike congruence x a 
is "comoving" ( "frozen-in" ) along the worldlines of the fundamental observers u a . In addition this fact 
ensures that T ab (x.) and 7?. a t,(x) lie in the screen space and the unit vector fields {x a , u a } are orthogonal 
at any spacetime instant. 

The incorporation of the EFE is achieved by applying (2.3) and (2.11) to the unit vector field x a . This 
leads to propagation (along x a ) and "constraint" (lying on the screen space) equations which describe 
the influence of the gravitational field in the spatial variation of the kincmatical and dynamical variables 
of the timelike and spacelike congruences. In the next section we shall derive the complete set of the 
dynamical equations associated with the orthonormal tetrad {u a , x a , y a , z a } in order to analyze and 
express the evolution and constraint equations (2.6)-(2.8) and (2.12)-(2.15) in terms of the kinematical 
quantities of the spacelike congruences with a view to reveal the special structure of the new constraint 
(2.15). 
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2.3 SIIS models and spacelike congruences 

From equations (2.20) and after appropriate projections, it follows that each of the eigenvectors {x a , y a , z a } 
has the following properties: 

x a = 7V a ( x ) = = ft ab (x) (2.30) 

Va = N a (y) = = TZ ab (y) (2.31) 

z a = N a (z) = = TZ ab (z) (2.32) 

which implies that each generator of the spacclikc congruence is comoving with the fundamental observers 
and the corresponding screen spaces are integrable submanifolds of M.. Projecting equation (2.8) or (2.15) 
with x a x b , y a y b and z a z b we obtain: 

T Qb (x) = a (y a y b - z a z b ) (2.33) 

Tab{y) = (i (z a z b - x a x b ) (2.34) 

T ab (z) = 7 (x a x b - y a y b ) (2.35) 

where a,f3, 7 are the eigenvalues of the associated shear tensors of the spacelike vector fields x a ,y a ,z a 
respectively. 

Furthermore projecting (2.8) and (2.15) along y a z b , z a x b and x a y b , the shear and electric part con- 
straints can be expressed as spatial variations of the corresponding eigenvalues along the individual 
spacelike curve: 



(03 - cr 2 )* = — ^ (cr 3 - o- 2 ) - 3aa 1 


(2.36) 


£ 

(ci - 0-3)' = — ± (ai - <7 3 ) ~ 3/3<t 2 


(2.37) 


(02 -Cl) = ^ (f72 - CTl) - 37(T3 


(2.38) 


(S 3 - £2)* = - y (#3 - E 2 ) - ZaE 1 


(2.39) 


(E 1 - E a )' = -^-(E 1 - E 3 ) - 3(3E 2 


(2.40) 


(E 2 - Eh) = -y (E 2 - E,) - 3 7 £ 3 


(2.41) 


where we have introduced the notations: 




= A„.. :i //. (K a ...) = K a .., k z k . 


(2.42) 


The advantage of decomposing the constraints along the specific spacclike curve is immediately apparent 
from equations (2.36)-(2.41): the associated spatial variations of the shear and electric part eigenvalues 
have been partially decoupled and expressed in terms of the corresponding kincmatical variables of the 
spacclikc congruences. Obviously, this decoupling process is the result of the conditions (2.30)-(2.32). 

The above set of equations must be augmented with the evolution equations (2.6), (2.12), (2.13) and 
the divergence equations (2.7), (2.14): 


a 


(2.43) 


E a = —-p& a — 0E a + 3a a E a — E a a a 

a 


(2.44) 


p = -pe. 


(2.45) 
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P=3 
p' = 3 
p = 3 



E' 2 



_3£x 
2 

3£ ±1 
2 



a(E 2 


-E 3 ) 


, 8= 


P(E 3 


-Si) 


, <?' = 


7 (Ei 


-E 2 ) 


, 9 = 



o"i — a (o-2 — 0-3) 



3£ 



3£. 



(72 -13 ((73 - CTl ) 



<7 3 + -^CT3 -7(°1 - CT 2) 



(2.46) 
(2.47) 
(2.48) 



Using a similar splitting of the Ricci identities for each spacelike vector field of the orthonormal triad 
{x a , y a , z a }, we obtain evolution (along u a ), propagation (along the individual spacelike vector field) and 
"constraint" (lying in the associated screen space) equations of the kinematical quantities of the spacelike 
congruences which we give them below. 

The evolution and the "constraint" equations arise from the temporal and the spatial projections of 
the trace part of the Ricci identity respectively: 



(£y) 
(SO" 



a (cr 2 - cr 3 ) , 



~ 7^£yO -0(<T3-O-l), 



£, 



a 



1^+13 



1 1 n 1 . 

a = -aai - -a8 - -£ x (ct 2 - cr 3 ) 



19=5^3-^-^(^8-^1) 

7 = ^70-3 - ^7# - (0-1 - 0-2) 



f 
~~2 



[|-7 



(2.49) 
(2.50) 
(2.51) 

(2.52) 



T~ 7 



( y + a 



[y + 7 



f +7) —*, *-/» 



(2.53) 
(2.54) 



It is interesting to note that the evolution equations (2.49)-(2.51) imply for example a = E x = for a 
Petrov type I model. In addition, by taking appropriate linear combinations of the "constraint" equations 
(2.52)-(2.54) this set is essentially equivalent with the Jacobi identities or the twist-free property 1Z = 
of the spacelikc congruences. 

The propagation equations of the expansion and the shear of the spacelikc curves follow from the trace 
and trace free part of the Ricci identity [44] which in the case of the SIIS models take the form: 



f 2 

-^-2a 2 
2 



f+7 



2p 



29 

+ 3)U 



(t+ 7 K(t 



0-1 
£ 



(2.55) 



(£y)' = 



£ 2 



~2 



Eo 



'ill 
3 



0-2 



a £ x 



26 



0-2 



i)£z 



(2.56) 



(£ z 



£ 2 



a3 + 3 



2p 



"7j °3 



— ex.] £ x 



(2.57) 



(a)* = 



F 1 ( E V 

-at x — 

2 I 2 



E x + 2E 2 1 . 



y + 7 ] 7 



(2.58) 



i?2 + 2_E.s 



7 + 



y-7 7+2 



- (cr 2 + 2(73) ( cr 2 , 



2 



2 



(2.59) 



(7 - - 7& -i(|- a )* + (|-«)« + i(| +/3 )' + 

-^-^ + -(.3+2.0 ^3 + 3 J 

The set of equations (2.36)-(2.41) and (2.43)-(2.60) completely characterizes the dynamics of SIIS models 
and it will be used in order to prove the conjecture by showing that vacuum and non- vacuum silent models 
of Petrov type I always admit a three dimensional group of isometries thus reducing to the subclass of 
SH models of the Bianchi type I. 

Finally it should be remarked that there exists a direct correspondence between the description of SIIS 
models with the theory of spacelike congruences and the Orthonormal Frame (ONF) approach [47, 49]. 
For example, the spatial curvature quantities ni2, n^x, "23, 0,1, Q-2, 0,3 and the kincmatical quantities of the 
spacelike congruences are related via: 

ni2 = -7 (2.61) 

«3 = -y (2.62) 

whereas the set (2.36)-(2.41) and (2.43)-(2.60) is the covariant version of the tetrad equations presented 
in [22]. As we shall see in the next section, this fact (already emphasized in [22]) makes clear the covariant 
property of the incompatibility of the constraints (2.36)-(2.41) with the evolution equations. At the same 
time reveals in a transparent way, the geometric nature of (2.36)-(2.41) since they induce, through the 
kincmatical variables of the spacelikc congruences, a set of symmetry constraints leading to the subclass 
of SH models of Bianchi type I. 



(2.60) 



"■23 



"■31 



-/3, 



a i = -y. a 2 = - — , 



3 Uniqueness of the Szekeres solutions within the family of SIIS 
models 

As we have mentioned in the Introduction, the vanishing of the magnetic part of the Weyl tensor gives rise 
to a set of constraints which in the 1+1+2 covariant formalism takes the form of equations (2.36)-(2.41). 
In order to be consistent with the evolution equations (2.43)- (2.45), they must not generate a new set 
of constraints after their repeated propagation along the timclike congruence u a . The outcome of this 
procedure has led to the proof of the consistency of Petrov type D models and that type I models develop 
a chain of new highly non-linear and non-trivial algebraic constraints that need to be solved [20, 22, 23]. 

However, it will be enlightening if we choose an alternative direction to deal with the consistency 
checking, in the spirit of the notion of geometric symmetries. The concept of a geometric symmetry 
or collineation admitted by the spacetime manifold, is closely related to that of a transformation that 
leaves a specific geometric object invariant (up to a conformal factor) along the integral lines of the 
generating vector field. Although there exists (and can be defined) a sufficiently large number of geometric 
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symmetries, they have, in general, the drawback of imposing severe restrictions on the geometry and 
the dynamics leading, in most cases, either to physically unsound models or models with very special 
properties (see for example [48] for SH models of Bianchi type I). In order to demonstrate this fact, let 
us consider the case of a spacelike Conformal Vector Field (CVF) X a = Xx a which is defined according 
to the relation: 

Cxg a b = 2ipg a b (3.1) 

where ip is the scale amplification or dilation of the spacetime metric. 

In terms of the kinematical variables of the spacelike congruence the conformal equation (3.1) can be 
shown to be equivalent to the conditions [46]: 

1 * 
T ah (x) = 0, {\nX). a = -E x x a - x a 

x a u a = ~S X , N a = -2tu ab x b (3.2) 

and the conformal factor is then given by ip = -^.XE X =X- 

Furthermore, the coupling of the geometry with the matter fields, through the EFE, implies that 
there is a mutual influence between any geometrical or dynamical constraint with subsequent restrictions 
in the structure of the corresponding models. For example, in the case of an irrotational and geodesic 
fluid flow the conditions (3.2) become: 

T ab (x)=0, (\nX). a = -x a 

S x = 0, N a = 0, V = 0=1 (3.3) 

i.e. the spacelike congruence is necessarily comoving with the observers u a and the CVF reduces to a 
KVF. Using the covariant description of the Petrov type D models in terms of the kinematical quantities 
of the spacelikc congruences (sec below), this in particular implies that the Szekeres solutions cannot 
admit spacelike CVF parallel to any of the eigenvectors of the shear (or the Weyl electric part) tensor. 

The consistency of the Petrov type D models along with their non-symmetry property (non-existence 
of isometries), allows one to speculate about the existence of some type of link between the consistency 
conditions and the presence of a geometric symmetry in the case of the Petrov type I models. This 
expectation can be justified by the fact that the new algebraic constraints in Petrov type I models are 
coming from the time propagation of the initial integrability conditions (2.39)-(2.41). In addition, they 
involve the kinematical and dynamical variables and necessarily 1 lead to a restriction of the underlying 
geometric structure of the corresponding models. On the other hand the existence of any geometric 
symmetry is controlled by a set of integrability conditions which, after repeated covariant differentiation, 
also lead to a system of algebraic relations. Consequently one may argue that the presence of the non- 
trivial constraints (2.39)-(2.41) (together with the rest of conditions imposed by (1.1)) in SIIS models of 
Petrov type I can be loosely interpreted as corresponding to some type of symmetry restriction. 2 In the 
following we shall prove that this is the case and the corresponding symmetry constraints arc equivalent 
to the Killing conditions (3.3). 

The above arguments and the conclusions they lead to, can be summarized in the following: 



1 Evcn if we assume that at most one, in the whole set of algebraic relations, is linearly independent this means that 
one of the kinematical or dynamical variables of the model can be expressed as a function of the other, thus leading to an 
associated restriction of their broadness. 

2 The interpretation of the induced constraints in terms of a set of symmetry conditions indicates a possible connection 
between the full set of EFE and the existence of a specific symmetry. For a subclass of SH vacuum and non-titled perfect fluid 
models, this connection has been established showing that a specific symmetry may play a role in the invariant description 
of general cosmological models [50]. 
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Theorem There do not exist vacuum and non-vacuum SIIS models of the Petrov type I with vanish- 
ing cosmological constant 



Before we proceed with the proof of the Theorem, it will be useful to covariantly characterize the SIIS 
models of the Petrov type D (So) in terms of the kinematical variables of the spacelike congruences. 
Assuming that E\ = E 2 a\ = a 2 , equation (2.38) (or (2.41)) implies that 7 = <^=> T a b(z) = 
p a k (z)p b l (z)zk ; i = ^fpab{%)- Therefore SIIS models of the Petrov type D are described by the conditions: 



Ei = E 2 0-1 = o- 2 , T ab (z) = 0, 



(3.4) 



It is interesting to note that under the conditions (3.4) the complete set of equations (2.36)-(2.41) and 
(2.43)-(2.51) is consistent under repeated propagation along the timelike congruence u a . Furthermore, us- 
ing the local coordinate form of the SIIS metric (2.21) and choosing x a = e^ A S^, y a = e~ B Sy, z a = e~ r <5" 
we can verify that the conditions (3.4) follow from the assumption A = B + H where the function H 
satisfies H = H = 0. 



Proof 

Consider the case of Petrov type I models (Si) which implies that E a ^ Ep and a a ^ op for every 
a, j3 = 1,2,3. Using the Ricci identity (2.3) it is straightforward to show the following commutation 
relations for every scalar quantity S (see also equation (A8) in [21]): 



S) = [S] -S[a 1 + - 



(3.5) 



(S') 



S) -S'{ (72 



- S [ 0-3 + - 



(3.6) 



(3.7) 



In order to analyzing the consistency of the constraints (2.39)-(2.41), we propagate them along u a and 
use the commutation relations (3.5)-(3.7) and the set of equations (2.43)-(2.44) and (2.46)-(2.51). As 
expected, the dynamical conditions (2.39)-(2.41) are not identically satisfied and they lead to the following 
new constraints: 



6 E2 (0-3 - 0-2) + 0-2 (E 3 - E 2 ) 



3£ x [E 2 (2a 2 -0-3)- E 3 a 2 ] 



-2a \E 2 (2cr 2 - II0-3) - E 3 ( llcr 2 + I60-3)] 



(3.8) 



6[E' 3 (a3-a 1 )+a 3 (E 3 -E 1 )} = 3S y [-E 3 (2a 3 - *i) + Ei<r 3 ] 



-2/3 [-Ea (2(t 3 - llcri) + E x (ll<r 3 + I601)] 



(3.9) 



6 E x {p<x-o-x)+ox(E 2 -Ex) = 3S Z [Ei (2a x - tr a ) - E 2 ax] 



-2j [Ex (2(7! - 11(72) - E 2 (11(7! + 16(7 2 )] 



(3.10) 



Clearly, equations (3.8)-(3.10) represent the covariant form of the constraints (67)-(69) given in [22]. 
Moreover, the assumption of a Petrov type D model (for example Ex = E2 <J\ =0-2) implies that 
7 = and the rest of equations are trivially satisfied. 
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Solving equations (3.8)-(3.10) with respect to E2, ^3? a second temporal propagation leads to the 
following expressions for the spatial derivatives cr 2 , a' 3 and ai: 



where: 



CT 2= 0- 3 = — , (Tl = — (3.11) 

A 2 £>2 <^2 



A 1 = -2E 2 2 [3£ x a 2 + 2a (a 2 - 10(7 3 )] 

+E 2 {4E 3 [3£ x a 2 + 2a (cr 2 + 8a 3 )] 

+3 (cr 3 - o- 2 ) [£ x (lOaf + 4(72(73 - 5of ) + 2a (l4of + 20cr 2 (7 3 + llcrf)] } 

-2J5$ [3£> 2 + 2a (19(72 + 8ct 3 )] 

-3E3 (£ x + 2a) ((7 2 - (7 3 ) (2of - 7(7 2 ct 3 - 4of) 

-/? (0-3 - cr 2 ) [3£ x a 2 (cr 2 - <r 3 ) + 2a (a 2 - 11(7 2 (73 - 8<r|)] 

A 2 = 6{2£; 2 2 -£;2[4-B3 + 3 ((73 -(7 2 )(2(72 +a3)] 

+2£ 3 2 + 3E 3 ((73 - <r a ) (2<r 3 + a 2 ) - p (a 2 - a 3 ) 2 }} 



Bx = -2^ 3 2 [3£ y (73 + 2/3((73-10(7 1 )] 

+E 3 {4Ex [3S y a 3 + 2(3 (<t 3 + 8<n)] 

+ 3 ((7l - (7 3 ) [£ y (10(7 2 + 4(7K7 3 - 5(7 2 ) + 2/3 (14(7 2 + 20(71(73 + ller 2 )] } 

-2E\ [3£ y a 3 + 2/3 (19(73 + 8ai)] 

-3Ej + 2/3) ((7 3 - (7i) (2(7 2 - 7(71(73 - 4(7?) 

-P (oi ~ 0-3) [3£ y (7 3 (a 3 - a,) + 2/3 ((7 3 2 - lltr^ - 8(7 2 )] 

B 2 = Q{2El-E 3 [AE 1 +i{a 1 -a 3 ){2a 3 + a l )] 

+2E'f + 2>E X (<n - (7 3 ) (2ai + o- 3 ) - P ((7i - (7 3 ) 2 ]} 



d = -2£ 2 [35,<7i + 2 7 (ffi - 10(7 2 )] 

+£i{4£ 2 [3f,o-i + 2 7 ((7! + 8(7 2 )] 

+3 ((7 2 - (7i) [£* (10(7 2 +4(7i(72 - 5(7 2 ) + 2 7 (14(7 2 + 20(7i(7 2 + H(7 2 )]} 

-2£ 2 [3£ z (7i + 2 7 (19(7i + 8(7 2 )] 

-3£ 2 (£, + 2 7 ) ((71 - (7 2 ) (2(7 2 - 7(71(72 - 4(7 2 ) 

-p ((7 2 - <T\) [3f 2 (7i ((7i - (7 2 ) + 27 ((7 2 - ll(7l(7 2 - 8(t|)] 

C 2 - 6{2£ 2 -£a[4£;2 + 3((72-(7i)(2(7i+(72)] 

+ 2£f + 3£ 2 (C7 2 - (71 ) (2(7 2 + (7l) - p ((71 - (7 2 ) 2 ]}. 

The above equations correspond to the covariant form of the 1+3 orthonormal frame expressions (72)-(73) 
of [22]. 
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In the reduction to the Petrov type D models i.e. as one takes the "limit" Si — > Sd by successively 
substituting 7 = and E 3 = — 2E 2 ,a 3 = —2a 2 , the new constraints must lead smoothly to identically 
satisfied relations. Therefore the consistency (<^> existence) of SIIS Petrov type I models implies that, in 
this "limit", the constraints (3.11) (and also equations (3.8)-(3.10)) will be reduced to the corresponding 
constraints of the Petrov type D models. Indeed we can verify that for Si — > Sd, equations (3.8)-(3.10) 
and (3.11) lead to the set (2.36)-(2.41) when we specialize it to Sd models. 

On the other hand the vanishing of the polynomials A, B,C in the equations (3.11) gives relations 
between the kinematical variables of the spacelike congruences of the form a ~ £ x , (3 ~ £ y , 7 ~ £ z . At 
the "limit" Si — > Sd these relations reduce to a = £ x /2,(3 = —£ y /2,^ = 0. We observe that the third 
equation reproduces the Petrov type D property 7 = 0. The rest of the "limiting" values suggest that in 
this case, if a Petrov type I model exists it should be reduced ( in that "limit ") to the subclass of Petrov 
type D models satisfying a — £ x /2 and (3 = —£ y /2. In fact we can easily show that this case corresponds 
to the Locally Rotationally Symmetric (LRS) models of Ellis class II [47, 51]. 

A further propagation of (3.11) and the use of (2.43)-(2.45), (2.49)-(2.51) and (3.5)-(3.7) give algebraic 
relations of the form: 

a = £ x f 1 (E 2l E 3 ,9,p,a 2 ,o- 3 ) (3.12) 



P = £yf2(Et,E 3 ,6,p, 0-1,0-3) 
7 = £zh(Ei,E 2 ,0,p, 0-1,02) 



(3.13) 
(3.14) 



where /i,/2,/3 are rational functions in which both the numerators and denominators are fourth or- 
der polynomials with respect to the factor E 2 — E± and their corresponding coefficients are fifth order 
polynomials in 02 — <J\ . 

As expected, the propagation of equations (3.11) is also not identically satisfied which means that 
(3.12)-(3.14) represent another set of non-trivial constraints. Proceeding in a similar way as before, we 
take the "limit" Si — > Sd and expect to obtain identities since, at this level of differentiation, there are no 
corresponding equations for the Petrov type D models. However, it can be shown that (with an obvious 
abuse of notation): 

lim fi = l, Jim f 2 = -\, „lim f 3 = 0. (3.15) 



Si—>Sn 



At this stage, the values of the "limits" (3.15) can be understood due to the specific structure of the 
rational functions (3.12)-(3.14). In particular, the linear parts of the polynomial expressions in f\ and f 2 
are proportional to each other with coefficients 1/2 and —1/2 respectively. The ratio of the associated 
linear parts in f 3 scales as ~ (£3 + 2E 2 )(a 3 + 2a 2 ) which implies that f 3 — > for Si — > Sd- Moreover, 
the indeterminate cases where either the two pairs of numerators and denominators or the linear parts of 
fi, f 2 vanish, lead to algebraic relations 8 = 9{a\, ...). However, it can be easily verified that none of these 
solutions is consistent with the evolution equations (2.5)-(2.6) and (2.12)-(2.13) except if A ^ 0, which in 
turn implies the existence of SIIS models of the Petrov type I with non-vanishing cosmological constant 
[42]. In the Appendix, we give the solutions of the indeterminate cases and for illustration purposes, we 
demonstrate the inconsistency of one of these solutions when the cosmological constant vanishes. 

Assuming that £ x ,£ y ,£ z ^ 0, from the evolution equations (2.49)-(2.51) and for every Petrov type I 
model, we obtain the following equations: 



(AT 



f ) ^(l-4/ 1 2 )(o 3 -o 2 )/4 



(l-4/f)(oi-o 3 )/4 



= (l-4/ 3 2 )(o 2 -oi)/4 



In 



In 



In 



2/1 


+ 1 


2/1 


- 1 


2/2 


+ 1 


2/2 


- 1 


2/3- 


f 1 


2/3 


- 1 



= o 3 — 02 



0-1 — 0-3 



02 — 01. 



(3.16) 

(3.17) 
(3.18) 
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It should be noted that in the last equations, we have calculated the time derivatives of the functions f a 
by using the evolution equations (2.49)-(2.51). In this sense, the initial constraints have been time prop- 
agated only up to third order, avoiding to compute the derivatives directly from the highly complicated 
expressions (3.12)-(3.14). 

Expressing equations (3.16)-(3.18) in (local) coordinate form, a straightforward calculation yields: 

e T+D 1 (x,y,z) ^ e B 
h = 2 (e r +- D i(*>^) ±e B ) 



e A+D 2 (x, v ,z) -p e r 

2( e A+D 2 (x,y,z) ± e r) 
B+D 3 (x,y,z) ^ £ A 



(3.20) 



2 ( e B+D 3 {x,y,z) ± e A) ^ 3 - 21 ) 

where D a (x,y, z) are smooth functions of the spatial coordinates. 

Equations (3.19)-(3.21) imply that for Sj -> S D => A = B + H we have / 3 = => D 3 {x,y,z) = 
H(x,y). However fx ^ 1/2 and f<x ^ — 1/2 for arty smooth functions A,B,Di,Ds. This contradiction 
shows that the a;, y expansion rates must vanish £ x = £ y = which by means of equations (2.49)-(2.50) 
implies that a = (3 = 0. Hence, due to the Killing conditions (3.3), there exists a (necessarily) Abelian 
group of isomctries acting on 2-dimcnsional spacelikc orbits. Furthermore, equations (2.30) and (2.31) 
imply that each of the spacelike KVFs is hypersurfacc orthogonal, therefore the SIIS models of Petrov 
type I are reduced to the subclass of diagonal G 2 models [39] . Using a similar procedure we can further 
show that £ z = 7 = 0. Consequently, from equations (3.3) or, equivalently, the correspondence (2.61)- 
(2.62), it follows that the Petrov type I model always admits three independent and commuting spacelikc 
KVFs and the resulting silent model being SH of the Bianchi type I. 

It should be remarked that the above methodology cannot be applied when: a) two 3 of the functions 
f a are equal to ±1/2 or b) when fx = f% = and fs = ±1/2 (a special class of the diagonal G 2 family of 
perfect models). Nonetheless, both cases are incompatible with the full set of the dynamical equations. In 
particular, in the diagonal G 2 class of models, satisfying fs = ±1/2 7 = ±£ z /2, equations (2.55)-(2.60) 
imply that 

9 2 - 30<7i-3(p + 3<7i<7 2 + 3<r|) 20 2 +30<7 2 -6/>-9<rf 
Ex = , E 2 = . (3.22) 

However, from (2.38), (2.41) and (2.48) it follows that the quantities w\ = cr 2 + 0/3 and w 2 — E 2 + p/6 
are spatially homogeneous. Using the commutation relation (3.7), we find that the algebraic equation 
(w 2 ) = has as a solution o~ 2 = a 2 (wi, w 2 ,ax). Substituting in the initial relation ^3 = ±1/2 we finally 
compute <7i = <7i{wx,w 2 ) and a 2 = o 2 (wx,w 2 ) i.e. the shear eigenvalues are also spatial homogeneous 
which in turn implies that £ z = = 7. Using a similar method we can easily show that the conditions 
fx = 1/2 and f 2 = —1/2 lead also to £ x = £ y = = a = j3 and should be excluded. 

We conclude this section by noticing that the above "limiting process" is always possible within the 
SIIS solutions of the EFE due to their specific geometrical and dynamical behaviour. In particular, 
let us assume that a Petrov type I model exists and is completely disconnected (in the spirit of the 
above reduction) from the Petrov type D models (for example, due to the occurrence of denominators 
of the form (E 2 — Ex) etc.). Then, this must be true at any point of the spacetime manifold or, strictly 
speaking, during the time evolution of the model. Nevertheless, it has been shown that (apart from a 
set of measure zero) for initial states within the Petrov type I invariant set, the orbits approach at the 
asymptotic regimes (i.e. at early or late times), either a Petrov type D model or a vacuum SH model of 
Bianchi type I [15, 16, 25]. The exceptional initial states are treated similarly. From equations (3.16)- 
(3.18) it turns out that f a arc monotone functions of the time coordinate t or the dimcnsionlcss time 
variable r which is defined according to dt/dr = 3/9. In addition and taking into account the relations 



3 Using the Ricci identities (2.52)-(2.54) for the spatial triad {x a ,y a , z a } it can be shown that e.g. a = E x /2 implies 
S x = = a or P = -£ y /2. 
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(3.16)-(3.18), the functions f a satisfy f a — > ±1/2 as r — > ±00 which due to equations (2.52)-(2.54) and 
the previous discussion implies that £ x = £ y — £ z = and the model will approach a SH model of Bianchi 
type I. 

4 Discussion 

In the present article we have undertaken a detailed covariant analysis of the question regarding the 
broadness of vacuum and non-vacuum SIIS models which is equivalent to the consistency of the spatial 
div and curl equations with the evolution equations. By exploiting the basic tools from the theory of 
the spacelike congruences and the fact that the shear a a b and the electric part E a b share a common 
spatial eigenframe, we have expressed the constraints (2.36)-(2.41) in terms of the spacelike irreducible 
kinematical variables. These expressions make clear the significant role that the spacelike expansions and 
shear magnitudes play in the sense that they influence the overall spatial variation of the kinematical 
and dynamical variables (and vice versa) of the corresponding models. Therefore in order to have a 
transparent and complete picture of the dynamics of SIIS models, we have incorporated the associated 
time and spatial propagation equations (2.49)-(2.60) of the spacelike kinematical quantities. 

Obviously, the suggested procedure has strong similarities with the 1+3 ONF approach and the 
basic steps in analyzing the SIIS are essentially the same, which shows the covariant nature of the 
incompatibility of the constraints (2.36)-(2.41) with the evolution equations. However, the approach 
followed here using elements from the theory of spacelike congruences, has allowed us to interpret the 
implications of the constraints (2.36)-(2.41) in a geometrical manner by expressing them as a set of 
conditions that must be satisfied by the kinematical quantities of the family of the spacelike curves. 

It is worth noticing that, since the spacclike vector fields have been chosen in a unique way (as 
eigenvectors of the shear tensor and the electric part of the Weyl tensor), these conditions represent 
invariant relations. This fact facilitated the analytical study of the induced constraints in a (local) 
coordinate system and allowed us to show that the SH models of Bianchi type I are the only compatible 
Petrov type I spacetimes within the family of silent models. 

We note that, the present approach permits us to draw further conclusions for more general spatially 
inhomogeneous configurations of the Petrov type I. For example, even when we allow the presence of 
a fluid with non-zero pressure (losing in that way the "silence" property), Petrov type D models are 
still the only permissible. This can be justified by observing that the pressure is necessary spatially 
homogeneous [5] (because x a = y a = z a = 0), hence when we apply the commutation relations (3.5)-(3.7) 
in the evolution of the eigenvalues of the electric part E a t (into which the fluid pressure enters) the form 
of the new constraints (3.8)-(3.10) is not affected. In addition, the fact that we didn't use equation (2.5), 
shows that the uniqueness of Petrov type D models is also independent, in general, from the presence of 
a non-zero cosmological constant. Nonetheless, an exception to this conclusion is when the eigenvalues 
of the shear tensor a a b are constrained to be proportional to 9, due to the appearance of indeterminate 
expressions. In this case, a non-zero pressure or the cosmological constant play a key role in the existence 
of a broader class of SIIS solutions apart from the SH models of Bianchi type I [42]. Finally we expect 
that this approach can be used in order to check how large is the family of SIIS models by relaxing 
the vanishing of the vorticity and the magnetic part of the Weyl tensor but maintaining their "silence" 
properties i.e. the vanishing of the curls of E a b, H a b and the pressure p. 

Appendix 

In this Appendix, we present the solutions of the indeterminate cases i.e. when either the two pairs of 
numerators and denominators or the linear parts of /i,/2 vanish. In particular, solving the algebraic 
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equations num(/i.2)= and denom(/i.2) = wc obtain 4 : 

(Ji (4Ei + 5E 2 ) 3ai (Si - £ 2 ) 



C7 2 



5J5i + 4£ 2 S^i + 4S 2 

2 [25gi 2 + Eh (4Qg 2 + 27(7?) + E 2 (l6E 2 + 27(7?) 
9af 



(4.1) 



a 1 (2E 1 + E 2 ) e _ 3ai (2#i+£ 2 ) _ 2 [£ff - 8£i£ 2 + £ 2 (lggg - 21af t 
4E 2 ~ E Y ' ~~ E 1 -^E 2 ' P ~ 9(7? 



<7 1 (4£; 1 - J B 2 ) 2 [£? + £?! (4^ - 27a?) + 4£f] 

a2= £; 1 + 2E 2 ' e = - 3ai ' p = 9^ ' 

On the other hand, the vanishing of the linear parts of f\, f 2 give the following solution 



2E 2 (E 2 - 3cr|) , . 

6 = -3a 2 , p= 2 ( 4 - 4 ) 

^2 

As we have mentioned in Section 3, none of the above solutions are consistent with the evolution equations 
when A = 0. For illustration purposes, we show the inconsistency of the solution (4.2). Propagating the 
second relation of (4.2) and using the evolution equations (2.5) and (2.43)-(2.45) we get the algebraic 
equation: 

E\ - 8£i£2 + 16-Ef - 9Act? = 

with solutions E\ = 4E 2 ± Sy/Aai. Obviously, these solutions arc valid only for non-zero (and positive) 
cosmological constant. Similarly, we can easily show the inconsistency for the rest solutions for the case 
A = 0. 
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